Coherent particle production in collisions of relativistic nuclei * 
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Here we give the results of our study of features of dense groups, or spikes, of particles produced in Mg-Mg 
and C-Cu collisions at, respectively, 4.3 and 4.5 GeV/c/nucleon aimed to search for a coherent, Cerenkov-like, 
mechanism of hadroproduction. We investigate the distributions of spike centers and, for Mg-Mg interactions, the 
energy spectra of negatively charged particles in spikes. The spike-center distributions are obtained to exhibit the 
structure expected from coherent gluon-jet emission dynamics. This structure is similar in both cases considered, 
namely for all charged and negatively charged particles, and is also similar to that observed recently for all-charged- 
particle spikes in hadronic interactions. The energy distribution within spikes is found to have a significant peak 
over the inclusive background, while the inclusive spectrum shows exponential decrease with two characteristic 
values of average kinetic energy. The value of the peak energy and its width are in a good agreement with those 
expected for pions produced in a nuclear medium in the framework of the Cerenkov quantum approach. The peak 
energy obtained is consistent with the value of the cross-section maximum observed in coincidence nucleon-nucleus 
interaction experiments. 



1. Introduction 

A coherent component of particle production 
mechanism, Cerenkov-like radiation, in high en- 
ergy particle collisions has been introduced a long 
time ago. The idea of mesonic Q and scalar (pi- 
onic) radiation has recently been a subject of 
the systematic analysis of production of mesons 
in high-energy pion-nucleon scattering in nuclear 
medium in terms of classical quantum mechan- 
ics H . Characteristic signatures of the Cerenkov 
mechanism, such as the differential cross-sections 
and angle-energy correlations of produced parti- 
cles, have been predicted. 

Another approach of the Cerenkov-like radia- 
tion in strong interactions was suggested within 
the QCD based coherent gluon-jet emission model 
Q. In this model the pseudorapidity distribu- 
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tions of centers of particle dense groups, called 
spikes, are proposed to be investigated. The 
distributions of spike centers are predicted to 
have two peaks due to destructive interference for 
quarks of the same colour (pp collisions) or to be 
singly peaked due to constructive interference for 
quarks of different colour (e.g. pp interactions). 
Recent observations in hadronic || interactions 
are found to be in agreement with these predic- 
tions. 

The aim of this report is to present our recent 
results on investigations of spikes in relativistic 
nucleus-nucleus collisions ^-^]. We search for a 
coherent component of hadroproduction mecha- 
nism studying two types of distributions assigned 
to the two above descibed scenarios. Respec- 
tively, on the one hand, we analyse the distri- 
butions of centers of spikes in the frame of the 
gluon-jet emission model and, on the other hand, 
we investigate energy distribution within spikes 
to find the features predicted by the nuclear pio- 
nic Cerenkov-like radiation (NPICR) approach. 

To note is that spikes have been extensively 
investigated last years using a stochastic picture 
of particle production, namely an intcrmittency 
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phenomenon has been searched for and obtained 
in all types of collisions Q . In our studies (!(]] , we 
also found the intermittency effect leading to mul- 
tifractality and, then, to a suggestion of a possible 
non-thermal phase transition during the cascad- 
ing. The latter observation have been confirmed 
in different reactions pH| . 

2. Experimental details 

The results presented here are based on the 
experimental data obtained after processing the 
pictures from the 2m Streamer Chamber SKM- 
200 @ installed in a 0.8 T magnetic field. The 
chamber was irradiated at the Dubna JINR Syn- 
chrophasotron. A beam of magnesium nuclei with 
momentum 4.3 A GeV/c was used to collide with 
a magnesium target, while a 4.5 A GeV/c carbon 
beam interacted with a copper target. A cen- 
tral collision trigger was used to start the Cham- 
ber if there were no charged or neutral projectile 
fragments (momentum per nucleon required to be 
greater than 3 GeV/c) emitted in a forward cone 
of 2.4°. A more detailed description of the set- 
up design and data reduction procedure are given 
elsewhere jl2]-[l4j]. Systematic errors related to 
the trigger effects, low-energy pion and proton 
detection, the admixture of electrons, secondary 
interactions in the target nucleus etc. have been 
considered in detail earlier and the total contri- 
bution is estimated to not exceed 3% |l^,|l5). 

A total of 14218 Mg-Mg events and 663 ones 
of C-Cu collisions were found to meet the above 
centrality criterion. In the utilized Mg-Mg sam- 
ple only negative charged particles (mainly 7r~ 
mesons with a portion of some 1% kaons) have 
been studied, while in the C-Cu sample all 
charged aprticles have been considered. 

In the Mg-Mg events the average measurement 
error in momentum (s p /p) was about 1.5% and 
that in the production angle determination was 
(e$) ~ 0.1°. The particles were selected in the 
pseudorapidity (r) = — In tan ^d) window of 
A?7 = 0.4 — 2.4 (in the laboratory frame), in which 
the angular measurement accuracy does not ex- 
ceed 0.01 in r\ units. The mean multiplicity of the 
selected pions in the Mg-Mg sample is 6.70±0.02. 

In the C-Cu sample the average measurement 



error in momentum (e p /p) was about 12%, the 
error in the polar angle was (e$) ~ 2°, and 
the angular measurement accuracy does not ex- 
ceed 0.1 in rj- units in the pseudorapidity range of 
A?7 = 0.2 — 2.6. In addition, particles with > 1 
GeV/c are excluded from the investigation as far 
as no negative charged particles were observed 
with such a transverse momentum. Under the 
assumption of an equal number of positive and 
negative pions, this cut was applied to eliminate 
the contribution of protons 

To overcome an influence of the shape of the 
pseudorapidity distribution on the results, we use 
the "cumulative variable" , 

V(f]) = r P(v')dv' / p(v')d V ' , (1) 

JVmin I ''imin 

with the uniform spectrum pW) within the inter- 
val [0,1], as advocated in Ref. Jlq] , This transfor- 
mation makes it possible to compare results from 
different experiments. 

The spikes are extracted in each event from 
the ordered pseudorapidities scanned with a fixed 
pseudorapidity window (bin) of size 5V- Spikes 
with a definite number of particles 5n, hit in the 
bin, are determined and, for each 5n, the distribu- 
tion of centers of spikes, averaged over all events, 
is obtained. The center of spike is defined by 

To reveal dynamical correlations, the Vq- 
distribution is compared with analogous distribu- 
tions obtained from the simulated pseudorapidity 
single-particle spectrum p^V) without any input 
information about particle correlations. The sim- 
ulation procedure was as follows. The number 
of particles was randomly generated according to 
the multiplicity distribution of the data sample. 
Then, the pseudorapidities were distributed in ac- 
cordance with the experimental ?7-spectrum cor- 
responding to the generated multiplicity. In each 
case of the reactions considered here, the total 
number of the simulated events exceeded the ex- 
perimental statistics by a factor of 100. It is clear 
that the statistical properties of this set are com- 
pletely analogous to those of an ensemble result- 
ing from arbitrary mixing of tracks from different 
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Figure 1. Experimental (•) and simulated (o) spike-center distributions in C-Cu collisions for different 
5?7-bins and multiplicities Sn: (a) Sv = 0.04, Sn = 4, (b) Sv = 0.08, Sn = 5, (c) Sv = 0.12, Sn = 7, (d) 
Sv = 0.2, Sn = 9. The curves represent Gaussian fits. 



events, subject to the condition of retention of 
the ^(^-distribution. So, the obtained sample 
represents independent particle emission. 

3. The results 

3.1. Spike-center distributions 

The pseudorapidity spike-center r] - 

distributions for four different-size <5??-bins and 
for spikes of various multiplicities Sn are shown 
in Figs. H and || for C-Cu and Mg-Mg collisions, 
respectively. 

For each reaction type, a two-peak structure of 
the measured distributions (solid circles) is seen 
with the peaks in the neighbourhood of the same 
TJo, independent of the width and multiplicity of 
spike. The shape of the distributions is in agree- 
ment with the structure predicted by the coherent 
gluon-jet emission model jl and is similar to that 
observed in hadronic interactions B . 



In order to estimate the position of the peaks 
and the distance between them, we fit these 
bumps with Gaussians and average over differ- 
ent spikes. The peaks are found to be placed at 
n ~ 0.17 and 0.57 corresponding to 770 = 0.60 ± 
0.05(stat) ± 0.12(syst) and 1.30 ± 0.03(stat) ± 

O.lO(syst) in C-Cu collisions, and in Vo ~ 0.19 
and 0.63 corresponding to 770 = 0.89±0.03(stat)± 
0.08(syst) and 1.63±0.05(stat)±0.10(syst) in Mg- 
Mg interactions. 

They are separated by the following do interval, 

d = 0.68 ± 0.06(stat) ± 0.16(syst) (C-Cu) 

d Q = 0.75 ± Q.06(stat) ± 0.13(syst) (Mg-Mg) (2) 

in r) units. These values are close to those from 
the above mentioned hadronic interactions. 
The dynamical origin of the structure obtained 

is seen from a comparison of the experimental ?7o- 
distributions with those based on the above de- 
scribed simulated events (open circles in Figs. @ 
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Figure 2. Experimental (•) and simulated (o) spike-center distributions in Mg-Mg collisions for different 
6^-bins and multiplicities Sn: (a) Sv = 0.05, Sn = 4, (b) SV =0.1, Sn = 5, (c) Sv = 0.15, Sn = 6, (d) 
Sv = 0.25, Sn — 7. The curves represent Gaussian fits. 



and g) . No peaks are seen in the simulated distri- 
butions, levelling off at the background, far below 
the measured peaks. This points to a dynamical 
effect in the formation of spikes in agreement with 
the coherent gluon radiation picture. 

In order to reveal a possible contribution of 
hadronic jets to the effect observed, for negative 
pions from Mg-Mg collisions we carried out ad- 
ditional analysis taking into account azimuthal 
collinearity. To this end, the alignment coeffi- 
cients, 



£IW cos 2(^-^0 

V n ( n - 1) 



(3) 



for an individual event with multiplicity n were 
used with being the ith particle azimuthal an- 
gle@. 

A [3 < criterion were applied to mini- 
mize a contribution of the jet-structural events. 

This reduced statistics twice. The resulted Vo- 
distributions are shown in Fig. || at the same Sv 



and Sn as in Fig. g without the /3-criterion. One 
can see that the structure of the ^-distributions 
for such azimuthally isotropic events does not 
change compared to that of Fig. [2] and the dis- 
tributions demonstrate two peaks. The Gaussian 
fits of the peaks give the peaks positions to be 
Vq « 0.88 and 1.63 with the distance between 
them d = 0.75 ± 0.06(stat) ± 0.15(syst). The ob- 
tained values are almost the same as those found 
from Fig. 0. 

The nearness of the positions of the peaks with 
and without /3-criterion points at the azimuthal 
isotropy of events with spikes. 

In order to assess the reliability of the results 
obtained, we studied the influence of the Ar- 
range used and the experimental error (e^) in the 
measurements of the polar angle Varying the 
A77 range and the error (e$) we found the struc- 
ture of the distributions unchanged and the po- 
sitions of the two peaks and the distance do to 
be within the above shown systematic errors, in 
support of the conclusions made. 
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Figure 3. The same as in Fig. ^ but for azimuthally isotropic events (see text and Eq. (^)). 



3.2. In-spike energy spectra 

The strong signal of the coherent emission dy- 
namics obtained allows further search for its man- 
ifestation in the energy distribution, as predicted 
in the NPICR approach ||. In this model, the en- 
ergy spectrum of pions, emitted through the co- 
herent Cerenkov-like mechanism when a few-GeV 
proton passes the nuclear medium, is predicted to 
have a peak. This peak is expected to appear at 
260 MeV when an absorption effect is neglected 
and at 244 MeV otherwise. 

In Fig. ^ we compare the c.m.s. in- 
clusive kinetic energy distribution, F(K*~) = 
(l/E*p*) dN/dK*, with analogous spectra calcu- 
lated for pions from spikes in Mg-Mg interactions. 
Here, E* and p* denote, respectively, the particle 
energy and momentum in the c.m.s. frame. 

Using the temperature description, utilized to 
characterize a system of excited hadrons |l4 , 18 
2 0| , we parametrize the inclusive spectrum (Fig. 
l|a) by a sum of two exponents, 

F(K*) = A, exp(-K*/T 1 )+A 2 exp(-iT/T 2 ) , (4) 

where the temperatures T\ < characterize [|l8| 
the two possible mechanisms of pion production, 



via A-resonance decay and directly, and are re- 
lated to the pion average kinetic energies. The 
range of the parametrization shown is limited 
from below and from above due to detector effects 
and corresponding requirements on the momenta 
of pions. The fit gives T x = 65 ± 1 MeV and 
T2 = 127 ± 1 MeV. These values are, in general, 
consistent with those obtained from the earlier 
analysis of the reaction under study |]l4]| and from 
other experiments (li|[nj . Some difference in the 
values could be explained if one takes into account 
the difference in the sizes of the (pseudo)rapidity 
regions used p0| . 

The shape of the F(K*) distribution changes 
when the analysis is extended to spikes, Figs. [|b 
and|c. The energy spectrum of particles belong- 
ing to a spike differs significantly from the expo- 
nential law (|J) and has a peaked shape. To ex- 
tract the NPICR-signal we compare the in-spike 
energy spectra with the renormalizcd inclusive 
distribution, or inclusive background, depicted by 
the dashed lines. The first peak is seen to be 
above background with a statistical significance 
of 2.7 and 4.1 standard deviations in Figs. [|b 
and 0c, respectively. This peak is located at the 
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Figure 4. Inclusive kinetic energy distribution (a) and analogous distributions for spikes of negative pions 
in Mg-Mg interactions, (b) 5V — 0.1, Sn — 6 and (c) Sv — 0.15, Sn = 7. The solid line represents the 



exponential fit with Eq. 14), the dashed lines show the inclusive background. 



kinetic energy K* w 100 MeV, or the total energy 
E* w 240 MeV, in accordance with the NPICR 
prediction. 

To estimate the position of the peak and its 
width and to make the results more compa- 
rable with the NPICR expectations, the E*- 
distributions of particles in spikes of various size 
(5??-bins and different <5n-multiplicities are stud- 
ied. Fig. H represents examples for these distri- 
butions. The following specific peculiarities are 
found. 

All these _E*-distributions possess a non- 
exponential behaviour with a pronounced max- 
imum in the vicinity of the value E^ = 240 MeV 
regardless of bin size and multiplicity of spikes. 
The higher the multiplicity of spike is (at fixed 

Sv size), the more peaks appear. A multi-peak 
structure is observed for bins with the multiplici- 
ties Sn > 3, while at Sn < 3 only one peak occurs 
(not shown). 

To reveal the dynamical signal we compare 



the in-spike energy distributions with the inclu- 
sive background (dashed-line). As for the above 
kinetic-energy distributions, the value of E* of 
about 240 MeV is obtained to be the position of 
the most prominently and statistical-significantly 
peak over background. To estimate the back- 
ground and to parametrize the signal, we use a 
fifth-order polynomial for the background and a 
Gaussian for the peak. The solid curve shows 
the result of this fit. After averaging over various 
spikes, the position of the peak and its width are 
found to have the values, 

E* m = 238 ± 3(stat) ± 8(syst) MeV, 

r m = 10 ± 3(stat) ± 5(syst) MeV, (5) 

respectively. 

The location of the obtained centre of the 
Gaussian lies within the interval of energies ex- 
pected for pions from the Cerenkov-like mecha- 
nism for incident protons of a few GeV, 224 < 
E m < 244 MeV g. The value of (§) is simi- 
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Figure 5. Total energy spectra for negative particles spikes of different SV-bms and multiplicities Sn 

in Mg-Mg collisions: (a) Sv = 0.05, Sn = 4, (b) 5rj = 0.08, Sn = 5, (c) Sv = 0.1, Sn = 5, (d) 

Sv = 0.15, Sn = 6. The solid lines represent the fit (see text), the dashed lines show the inclusive 
background. 



lar to the position of the peak observed in the 7r + p 
invariant mass distribution in the analysis of co- 
incidence measurements of (p,n) reactions on car- 
bon at 1.5 GeV/c in the A-resonance excitation 
region pi]] , the effect connected with the NPICR 
mechanism EJ. Also, the width T m confirms an 
observation of the Cerenkov radiation signal ex- 
pected to be T < 25 MeV. 

4. Conclusions 

In summary, in order to search for a coher- 
ent, Cerenkov-like emission mechanism of par- 
ticle production, a study of spikes in relativis- 
tic nuclear collisions is carried out with charged 
particles from central C-Cu collisions at a mo- 
mentum of 4.5 A Gev/c and with negative pions 
from central Mg-Mg collisions at a momentum of 
4.3 GeV/c per incident nucleon. The spike-center 
distributions and the energy spectra of particles 
within a spike are investigated for various narrow 
pseudorapidity bins and different spike multiplic- 
ities. 



The spike-center distributions are found to pos- 
sess a double-peak shape that is in agreement 
with the structure expected from the coherent 
gluon radiation model. The obtained distance 
between the peaks as well as the shape of the 
distributions are similar to those observed re- 
cently in analogous studies of charged-particle 
spikes in hadronic interactions. The dynamical 
effect in the spike-center distributions is revealed 
in a comparison with an independent particle- 
emission model, where no peaks are seen. 

The coherent character of the particle- 
production mechanism is confirmed by studying 
energy distributions in Mg-Mg interactions. The 
inclusive energy spectra show monotonic expo- 
nential decrease with two specific temperatures, 
while the in-spike energy distributions are ob- 
tained to exhibit a peak at a position and of a 
width both consistent with the values expected 
from the theoretical calculations based on the 
hypothesis of nuclear pionic Cerenkov radiation. 
The value of the peak energy is close to the 
recently observed maximum in the differential 
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cross-sections studied in the coincidence experi- 
ments of a few GeV nucleon-nucleus interactions. 

The results of the presented analysis signal- 
ize coherent emission as a complementary mech- 
anism to the stochastic scenario of hadroproduc- 
tion. Furthermore, the similarity of the spike- 
center distributions obtained for like-charged par- 
ticles in the presented paper with those for all- 
charged-particle spikes in the earlier studies in- 
dicates important contributions of the coherent 
mechanism to the formation of Bose-Einstein cor- 
relations p^ |. It is worth to mention that, in 
comparison to stochastic (intermittency) dynam- 
ics, the origin of which remains still unclear 0, 
the coherent emission has definite underlying dy- 
namics. All this gives evidence for the necessity of 
further efforts in studying existing experimental 
data. 
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